Uncollapsing the wavefunction by undoing quantum measurements 
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We review and expand on recent advances in theory and experiments concerning the problem 
of wavefunction uncoUapse: Given an unknown state that has been disturbed by a generalized 
measurement, restore the state to its initial configuration. We describe how this is probabilistically 
possible with a subsequent measurement that involves erasing the information extracted about the 
state in the first measurement. The general theory of abstract measurements is discussed, focusing 
on quantum information aspects of the problem, in addition to investigating a variety of specific 
physical situations and explicit measurement strategies. Several systems are considered in detail: 
the quantum double dot charge qubit measured by a quantum point contact (with and without 
Hamiltonian dynamics), the superconducting phase qubit monitored by a SQUID detector, and an 
arbitrary number of entangled charge qubits. Furthermore, uncoUapse strategies for the quantum dot 
electron spin qubit, and the optical polarization qubit are also reviewed. For each of these systems 
the physics of the continuous measurement process, the strategy required to ideally uncoUapse the 
wavefunction, as well as the statistical features associated with the measurement is discussed. We 
also summarize the recent experimental realization of two of these systems, the phase qubit and the 
polarization qubit. 
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I. INTRODUCTION 

The irreversibility of quantum measurement is an ax- 
iomatic property of textbook quantum mechanicsii In his 
famous article Law without law, John Wheeler expresses 
the idea with poetic flare: "We are dealing with [a quan- 
tum] event that makes itself known by an irreversible 
act of amplification, by an indelible record, an act of 
registration."^ However, it has been gradually recognized 
that the textbook treatment of an instantaneous wave- 
function collapse is really a very special case of what is in 
general a dynamical process - continuous quantum mea- 
surement. Continuous measurements do not project the 
system immediately into an eigenstate of the observable, 
but describe a process whereby the collapse happens over 
a period of timci^ The fact that continuous measurement 
is a dynamical process with projective measurement as a 
special case, leads us to ask whether the irreversibility of 
quantum measurement is also a special case. The purpose 
of this paper is to review and expand on recent develop- 
ments in this area of research, showing that it is possible 
to undo a quantum measurement, thereby uncollapsing 
the wavefunction, and to describe this physics in detail 
for both the abstract and concrete physical realizations. 



This paper follows our earlier work on the subject 
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as well as other papers investigating similar questions<^ 
Wavefunction uncoUapse teaches us several things about 
the fundamentals of quantum mechanics. First, there is 
a notion that wavefunction represents many possibilities, 
but that reality is created by measurement. The fact 
that the effects of measurement can be undone suggests 
that this idea is flawed, or at least too simplistic. If you 
create reality with quantum measurements, does undoing 
them erase the reality you created? Secondly, there is a 
wide-spread belief that quantum measurement is nothing 
more than a decoherence process. This suggests that the 



superposition never really collapses; it only appears to 
collapse. What actually happens, according to this idea, 
is that all the information about the system disperses 
into the environment: when a quantum system interacts 
with a classical measuring device, it becomes irreversibly 
entangled with all the particles that make up the mea- 
suring device and its surroundings. The uncoUapse of 
the wavefunction demonstrates that decoherence theory 
cannot be the whole story, because a true decoherence 
process is irreversiblei^ The perspective we take in this 
paper further advances the 'quantum Bayesian' point of 
view, where the quantum state is nothing more than a 
reflection of our information about the system. When 
we receive more information about the system, the state 
changes or collapses not because of any mysterious forces, 
but simply as a result of Bayesian updating. 

Our work is indirectly related to the 'quantum eraser' 
of Scully and Driihli^ There, the which-path information 
of a particle is encoded in the quantum state of an atom, 
resulting in a destruction of interference fringe visibility. 
On the other hand, if the which-path information of the 
particle is erased, the interference fringes are restored. 
Both the Scully proposal and our proposal erase informa- 
tion. However, there is an important difference. In order 
for the uncollapsing procedure to work, we have to erase 
the information that was already extracted classically. In 
the 'quantum eraser', only potentially extractable infor- 
mation is erased. 

While the first part of this paper deals with the ab- 
stract idea of uncoUapse, and formalizes its properties 
in terms of generalized quantum measurements, a great 
deal of the paper deals with actual measurement pro- 
cesses in specific physical contexts. This brings to mind 
the saying of Asher Peres: "Quantum phenomena do not 
occur in a Hilbert space. They occur in a laboratory" 
Following Peres' dictum, we discuss measurement pro- 
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cesses in a variety of solid state systems, where there has 
been remarkable experimental progress in recent years. 
Quantum coherence has been demonstrated to occur in 
a controllable fashion in systems such as semiconductor 
quantum dots and superconducting Josephson junctions. 
We will discuss the physics of measurement in these sys- 
tems, as well as concrete strategies for uncoUapsing the 
wavefunction. It should be stressed that two of these pro- 
posals (a superconducting phase qubit and optical polar- 
ization qubit) have now been implemented in the labora- 
tory, providing conclusive demonstration of wavefunction 
uncoUapse.^^ 

The paper is organized as follows. In Sec. II we de- 
scribe in detail what we mean by the undoing of a quan- 
tum measurement (Sec. II A), and give a general treat- 
ment of the physics, using the formalism of positive 
operator- valued measures (POVMs). This is done both 
for pure states (Sec. II B) and mixed states (Sec. II C). 
In Sec. II D we discuss the interpretation of wavefunction 
uncollapse, and what it tells us about quantum informa- 
tion. Sec. Ill begins discussion of the physical imple- 
mentation of this physics, with a treatment of the double 
quantum dot qubit, monitored by a quantum point con- 
tact. The measurement dynamics is discussed in IIII Al 
and we discuss the uncoUapsing strategy and statistical 
predictions in Sec. IIIIBI for the case of a qubit under- 
going measurement dynamics only. The results are dis- 
cussed in IIII CI and compared with the abstract results. 
In Sec. IIII Dj we discuss the statistical features of the time 
required to wait until the wavefunction is uncoUapsed. 
We generalize to the case of the finite-Hamiltonian qubit 
undergoing the uncollapse process in Sec. IIII El and de- 
rive results for the success probability in that case. We 
switch to a new physical system in Sec. IV, the super- 
conducting phase qubit. In Sec. IV A we discuss how 
the measurement process works for the phase qubit, and 
in Sec. IV B we discuss the uncoUapsing strategy, corre- 
sponding success probability, and the recent experimen- 
tal realization, which demonstrated the uncoUapsing. In 
Sec. V, we generalize to the case of many charge qubits, 
and discuss an explicit procedure to undo any general- 
ized measurement. We discuss recent developments in 
the theory and experiments of measurement reversal in 
Sec. VI and conclude in Sec. VII. 



II. GENERAL THEORY OF UNCOLLAPSING 
A. Preliminary discussion 

Our goal is to restore an initial quantum state dis- 
turbed by measurement. However, it is important to 
discuss what exactly we mean by that. For example, 
if we start with a known pure state \ipin) and perform 
a textbook projective measurement, then it is trivial to 
restore the initial state: since we also know the post- 
measurement wavefunction IV'm), we just need to apply 
a unitary operation which transfers \ipm) into \ipin)- If we 



start with a known mixed state, then its restoration af- 
ter a projective measurement is a little more involved)^ 
however such a procedure still can be easily analyzed 
using standard quantum mechanics and classical proba- 
bility theory. 

In this paper we consider a different, non-trivial sit- 
uation: we assume that an arbitrary initial state is un- 
known to us, and we still want to restore it after the mea- 
surement disturbance. To make this idea more precise, 
we consider a contest between the uncollapse proponent 
Plato, and an uncollapse skeptic Socrates. Socrates pre- 
pares a quantum system in any state he likes, but it is 
unknown to Plato. Socrates sends the state to Plato, 
who makes some measurement on the system, verified by 
the arbiter Aristotle. Plato then tries to undo the mea- 
surement. If Plato judges that the attempt succeeded, 
the system is returned to Socrates, with the claim that 
it is in the original state. Socrates is then allowed to try 
and find a contradiction in any way he likes, with the 
whole process monitored by Aristotle. If a contradiction 
can be found, then he can claim to refute the uncollapse 
claim, but in the absence of contradiction, the uncollapse 
claim stands. If Socrates would like to try again to find 
a contradiction, or if Plato judges that his undoing at- 
tempt was unsuccessful (and does not return a state), 
then Socrates prepares a new (still unknown to Plato) 
state, and the competition continues. If Socrates cannot 
find a contradiction after many rounds of the competi- 
tion, then Plato will win the contest, and will have suc- 
cessfully demonstrated the uncoUapsing of the quantum 
state. 

A slightly different but equivalent situation is when 
we know the initial state, but our uncoUapsing proce- 
dure must be independent of the initial state (so we can 
pretend that it is unknown), and therefore the uncol- 
lapsing should restore any initial state in the same way. 
This formulation is most appropriate for a real experi- 
ment demonstrating the uncoUapsing. Finally, we may 
consider the more general case where the measured sys- 
tem is entangled with another system, and we wish to 
restore the initial state of the compound system without 
any access to its second part. 

The traditional statement of irreversibility of a quan- 
tum measurement can be traced to the fact that it may be 
described as a mathematical projection. Projection is a 
many-to-one mapping in the Hilbert space, and therefore 
the same post-measurement state generally corresponds 
to (infinitely) many initial states It is therefore impos- 
sible to undo a projective measurement. 

However, the situation is different for a generali^ 
(POVM-type) measurement, which typically corresponds 
to a one-to-one mapping \ipm) — > IV"™) in the Hilbert 
space of wavefunctions (in this paper we consider only 
"ideal" measurements which do not introduce extra de- 
coherence). In this case the post-measurement wavefunc- 
tion IV'm) can still be associated with the unique initial 
state \ipin)-, and a well-defined inverse mapping exists 
mathematically. This makes the uncoUapsing possible 
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in principle. Since the inverse mapping is typically non- 
unitary, it cannot be realized as an evolution with a suit- 
able Hamiltonian. However, it can be realized using an- 
other POVM-type measurement with a specific ("lucky") 
result. 



B. Formalism for wavefunctions 

Let us first consider a pure initial state |'(/'m), and 
postpone a generalization to mixed states until Sec. 
lie. In the formalism of a general (ideal) quantum 
measurement^ which transfers pure states into pure 
states, the measurement with result m is associated with 
the linear Kraus operator M^, so that the probability of 
result m is 



(1) 



where 1 1 . . . 1 1 denotes the norm of the state, and the (con- 
ditioned) state after measurement is 



\lpr, 



v/^m(|^»„»' 



(2) 



where the denominator makes |?Am) properly normal- 
ized. (Very often people prefer to omit this denomina- 
tor and work with non-normalized states; this makes the 
mapping hnear.) The operators Em = M^Mm (called 
POVM elementsi^) are Hermitian and positive semidefi- 
nite by construction; these operators must obey the com- 
pleteness relation J2m^"i = 1, which ensures that the 
total probability of all measurement results is unity. A 
measurement operator Mm can always be written as 



M„ 



(3) 



where Um is a unitary operator (an important special 
case is when = ^/Em', this corresponds to the "quan- 
tum Bayes theorem"—). 

Now let us discuss wavefunction uncoUapse in this gen- 
eral and abstract context. The state disturbance rule ^ 
is typically a nonunitary one-to-one map in the Hilbert 
space. To undo the measurement with known result to, 
we have to realize a physical process corresponding to 
the nonunitary inverse operator M^^, multiplied by an 
arbitrary constant (which is not important because of the 
normalization). This can be accomplished with another 
measurement, possibly together with unitary operations. 
As shown below, we can realize measurement undoing if 
the second measurement realizes a Krauss operator of the 
form 



L ^ CULEm'^'^VL, 



(4) 



where Ul and Vl are any unitary operators, and C is an 
unspecified constant that will be discussed later. (The 
operator L also has a decomposition of the form UVe, 
but with a different POVM element E.) The uncoUapse 



then consists of three steps: first, the unitary operator 
Vlul^ is applied to reverse the unitary part of and 
prepare for the second measurement. Next the measure- 
ment operator L is applied. Finally the unitary operator 
is applied to reverse the remaining unitary part of L. 
We can now see the effect of the uncoUapsing operation 
on the state \ipm) by applying Eq. ([2]) and the unitaries 
to find 



\ulLvlul\4,m)\\ 



(5) 



thus restoring 
UlLVlulM, 



the original state, because 
m-^y-^m — C, which is removed by the nor- 
malization. (The phase of C is not important, since it 
affects only the overall phase of the wavefunction.) 

However, in order for the operator L to be physically 
realizable, the operator L must belong to another com- 
plete set of POVM elements, and therefore all its eigen- 
values must not exceed unity (otherwise some states will 
be assigned probabilities that are above unity; notice that 
the eigenvalues are non- negative automatically). Since 
L^'L = jCpy^-^m^^L, its eigenvalues are directly related 
to the eigenvalues p^™'' of the operator Em- Expressing 

Em — where the eigenvectors |z) form an 

orthonormal basis, the eigenvectors of L'^L are obviously 
Vl\i), and the corresponding eigenvalues are |Cp/p-'"\ 
Since all these eigenvalues must not exceed 1, we find the 
following inequality on |Cp, 



|Cp < minpj'™-' = minP,„, 

i 



(6) 



where min Pm is the probability of the result m, mini- 
mized over all possible states li^in) in the Hilbert space. 
The equality of min Pm to min^ Pj-"*^ follows from Eq. ([1]). 

There is no guarantee that the uncoUapse can be ac- 
complished deterministically, since we rely on a measure- 
ment with a specific result, corresponding to the operator 
L. We can calculate the uncoUapse success probability 
Ps from Eq. ([T]), with Mm L and lipm 



Ps = \\LVlul\ij, 



\c\' 



(7) 

Now using the bound ^ for |Cp, we find the bound 
for the success probability of uncoUapsing after the first 
measurement with result to: 



Ps < 



PJIi^'' 



(8) 



where the denominator is the probability of the result 
TO for a given initial state, while the numerator is this 
probability minimized over all possible initial states. 

The bound ([5]) is one of the most important results 
(notice a similar result in Ref. Q) and deserves a dis- 
cussion. First, this bound is exact in the sense that it is 
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achievable by the optimal uncoUapsing procedure. This 
is because the uncoUapsing operator with |C| = \J min Pm. 
is still a physically allowed operator. As we will see 
later, the upper bound ([5]) is achievable in real experi- 
mental setups. However, non-optimal uncoUapsing pro- 
cedures, especially involving a sequence of measurements, 
can lead to smaller success probabilities (an example 
of non-optimal uncoUapsing has been discussed in Ref. 
[l6j). An analysis of the procedures with an arbitrary 
sequence of measurements and unitary operations is sim- 
ilar to the above: the corresponding measurement and 
unitary operators should simply be multiplied. 

Notice that the success probability ([7]) and the inequal- 
ity ([5]) depends on the initial state, which is unknown to 
the person performing the uncoUapsing (Plato, see de- 
scription in the previous subsection). Therefore, the suc- 
cess probability Pg can be calculated by the man who 
knows what the initial state is (Socrates), while Plato 
can only estimate Pg; for example, he can calculate the 
worst-case scenario (the minimum of Pg over the accessi- 
ble Hubert space) or can calculate the average of Pg over 
all possible initial states (this procedure will be discussed 
in the next subsection). 

Recalling the fact that it is not possible to undo a fully 
collapsed state due to the nature of projective measure- 
ment, the uncoUapsing probability P5 should decrease 
with increasing strength of the first measurement. Qual- 
itatively, a stronger measurement is one that tends to 
a projection, as the uncertainty in the measurement de- 
creases. Mathematically, this means that some eigen- 
values of Em become closer to 0. As a consequence, 
minPm becomes smaller [see Eq. (|6])], therefore lowering 
the upper bound for P5. For a projective measurement 
Pg = minP„i = 0, thus making the state uncoUapsc im- 
possible. 

It is interesting to discuss the case when the initial 
state I V'in) is known to belong to a certain subspace of the 
Hilbert space, and we therefore wish to restore states only 
in this subspace. In this case, the calculation of min Pm 
should be limited to this subspace, which may increase 
the bound ([5]) for the success probabiUty P5. A trivial 
example of such a situation is when the initial state \tpin) 
is known to Plato. Here it is not necessary to minimize 
Pm over all possible initial states in Eq. ([5]), because the 
set of possible states consists of only one (known) state, 
thus allowing uncoUapsing with 100% probability. This 
is exactly the case discussed at the beginning of Sec. Ill Al 



C. Formalism for density matrices 

So far we have dealt only with pure states; however, it 
is very simple to generalize the uncoUapsing formalism to 
include density matrices. In this case the initial density 
matrix pi„ is transformed by the first measurement into 
the statei^ 



MmPiuMl^ 



where the probability P,„ of the measurement result m 
is 



Pm(Pm) = Tr(A'/^A'/„Pi„). 



(10) 



Using the uncoUapsing procedure previously discussed 
and using the same measurement operator L given by 
Eq. ([U, we find that the uncoUapsed state 



Pf = 



UlLVlulpmUmVlL^UL 

Tv{L^ LVLUlpmUmVl) 



= Pi: 



(11) 



coincides with the initial state. The uncoUapsing success 
probability Ps is equal to the denominator in Eq. (fTTjl . 
and satisfies the relation 



Ps 



\C\yPm{Pr. 



(12) 



as m Eq. dZl)]. The constant |Cp is stiU limited by the 
inequality ([5]), and therefore the probability of success 
has the upper bound 



Ps< 



iPrr. 



(13) 



(9) 



which is the same as the bound (O, except for the new 
notation in the denominator, which reminds us of the 
possibly mixed initial state. The minimization of Pm in 
the numerator should now be performed over the space 
of all possible initial mixed states; however, the result 
obviously coincides with the minimization over the pure 
states only. Similar to the case discussed in Sec. Ill Bl the 
inequality (fT3)) is the exact bound; it is achieved by the 
optimal uncoUapsing procedure, which maximizes \C\. 

If the initial state is pure, then the formalism of this 
subsection is trivially equivalent to the formalism of Sec. 
IIIBI It becomes more general in the case when the "ac- 
tual" initial state is mixed; for example, this happens 
when the initial state has been in contact with an un- 
monitored environment or Socrates prepares a state by 
a blind random choice from a set of pure states. A more 
interesting case for the result (jl3p is when the measured 
system is entangled with another system, which does not 
evolve by itself. Then the formalism can be applied to 
the compound system; however, the measurement proba- 
bility Pm depends only on the reduced initial density ma- 
trix, traced over the entangled second part. Therefore, in 
the entangled bipartite case the uncoUapsing procedure 
restores the state of the whole system, while the success 
probability P5 is given by Eq. ([13]) with pin being the 
reduced density matrix. 

Another advantage of Eq. (|13p in comparison with 
Eq. ([S]) is the following. In the derivation of both re- 
sults the initial state is the "actual" initial state, which 
is known to Socrates, but typically unknown to Plato. 
However, as we will prove below, Eq. (fT3)) can still be 
used by Plato in a somewhat different sense: with pin 
being understood as an averaged density matrix repre- 
senting a distribution of possible initial states. In this 
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case, Eq. p3)) gives the uncoUapsing probability aver- 
aged over this distribution. For example, if Plato knew 
that Socrates' strategy is to prepare one of 2 possible 
(nonorthogonal) states \tpi),\ip2), with probabilities V 
and 1 — V, then he could find the average undoing prob- 
ability in two ways. The first method is that he could 
simply average the undoing probabilities of the two states 
(also taking into account the information acquired in the 
first measurement, see below). Alternatively, he could re- 
call that the random state preparation described above 
is equivalent to considering the initial density matrix 
Prn = V\^i){^i\ + (1 - r)\^2){4'2\, and then applyil 
the result p3)) to this density matrix. In this way, in 
the absence of any information, Plato could estimate his 
typical success rate by calculating (fT3|) for a fully mixed 
state, invoking the principle of indifference 

In the general case the above statement, that both 
ways of computing the averaged undoing probability are 
equivalent, can be proven both logically and explicitly. 
For the logical proof we notice that Plato's judgment of 
successful undoing does not depend on whether or not 
Socrates knows the randomly picked state; therefore, the 
average probability of the cases judged to be success- 
ful should be the same in both situations (whether or 
not Socrates knows what the state is). Now let us also 
prove this statement explicitly, thus checking that our 
formalism is self-consistent. Suppose the initial state is 
prepared by Socrates by choosing randomly from a set of 
initial states p''^-' with probabilities Vk (the most natu- 
ral case is when initial states are pure, p'-'^^' = |V'fc)(V'fc|; 
however, this is not necessary). Then the bound for the 

average probability of uncoUapsing success P^""^ is the 
average of the bounds 



p(av) ^ \ ' 
- 



minP„t 



(14) 



Notice, however, that V^. is the posterior probability dis- 
tribution given the result m, which is different from Vk 
We may now invoke the classical Bayes rulei^ii^ 

Vim\k)Vk 



P(fc|m) 



Es^(m|fc)P^ 



(15) 



to relate the posterior V'/. = V{k\m) to the prior Vk and 
the conditional probability V{m\k) — Pm{p''^^) to have 
result m given state k, so that 



VI 



P,n{p^''^)Vk 
EkP^iP^'^)^ 



(16) 



Substituting Eq. (dH) into Eq. HH) and using J^k '^fe = 1 
in the numerator, we obtain 



p: 



(av) 



< 



min P„ 



minP„ 



(17) 



where = HkP^^^T^k is the averaged initial state. 

This ends the proof that Eq. (fTS]) can be used for an 
unknown initial state, with pj„ being understood as the 
average of all possible initial states. 



D. UncoUapsing probability, information, and an 
irreversibility measure 

We defined the success probability Ps as a probability 
to uncoUapse the post-measurement state pm- We now 
wish to start counting the overall success probability Ps 
from the time before the first measurement, so that Ps is 
the probability of the pair of events: measurement with 
result m and then successful uncoUapsing. Using Eqs. 
p2)) and p3|) we easily find the relation 



(18) 



Ps^PUp^n)Ps^\C? 

and the upper bound 

Ps < minPm. 



(19) 



Notice that Ps is independent of the initial state. While 
this property seems to be somewhat surprising, we will 
see later why it is rather obvious. 

If we now wish to consider all possible results of the 
first measurement, and perform different uncoUapsing 
procedures for each measurement result, then the total 
probability of uncoUapsing p^°*a-i is bounded as 



3total 



< 



min P„ 



(20) 



and is also independent of the initial state. 

The bounds (|19p and (|20p are exact and reachable by 
optimal uncoUapsing procedures. The bound (|20p indi- 
cates that 1 — minPm can be used as a measure of 
irreversibility (collapse strength) due to the measurement 
operation. 

Now let us discuss the relationship between the uncol- 
lapsing procedure and our knowledge of the initial state. 
While uncoUapsing is possible even if we know nothing 
about the initial state of the system, at first glance it 
seems like we gain some knowledge about the initial state 
in the process. This leads to the following interesting 
paradox, initially considered by Royer.^" By doing both 
a measurement and unmeasurement, one can seemingly 
learn something about the initial state without disturb- 
ing it. Then by repeating the measurement -I- unmeasure- 
ment process many times, even though the probability of 
such an event rapidly decreases to zero, the successful 
event would lead to essentially perfect knowledge of the 
initial state, leaving the state itself perfectly intact! One 
could then violate a host of known results, such as the 
no-cloning theorem. The resolution of the paradox lies 
in the fact that the pair of measurement and unmeasure- 
ment actually brings exactly zero information. Uncol- 
lapsing the state can only occur when the information 
in the second measurement exactly contradicts the infor- 
mation gained in the first measurement, thus nullifying 
it. This can happen in weak quantum measurements be- 
cause there is uncertainty about the system in the mea- 
surement result. It is to the extent that this ambiguity 
exists that it is possible to undo the weak measurement. 
Let us examine this in more detail. 
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We learn something about a pre-measurement state 
when the measurement resuh depends on the state. The 
measurement with resuh m brings some information 
about Pin because the probabiUty Pm{pin) depends on 
the initial state pi„. The ability to successfully uncol- 
lapse the state also brings some information about the 
initial state because the uncoUapsing probability Ps = 
\C\'^ I Pm{pin) also depends on pin- However, the collapse- 
uncoUapse probability Ps of observing both the result m 
followed by a successful uncoUapse is independent of pi^ 
- see Eq. ([T5| . Therefore, the combined effect of partial 
collapse and uncoUapse brings no information about the 
initial state. 

More quantitatively, we can use the same framework 
as at the end of Sec. 



II CI in order to track the informa- 
tion gain during the procedure. Suppose Plato assigns an 
initial distribution Vk of possible initial states as a statis- 
tical prior, to be updated as more information comes in. 
The measurement with result m brings in this informa- 
tion, so Plato updates his prior to the posterior distribu- 
tion V'f. [see Eq. P^ ]. Calculating in a similar way the 
distribution V'^ after the pair of the measurement and 
unmeasurement results, we find 



Ps{p^'^)Vk 



(21) 



where Ps{p) and Ps{p) denote respectively the prob- 
abilities of uncoUapsing (fT3)) and a combined coUapse- 
uncoUapse pair (fT5|) for the initial state p. However, as 
we have already stressed, Ps{p) is independent of the 
initial state p, and therefore cancels out of the expres- 
sion ((2T|) . This fact (and the normalization of the prior 
{Pk}) restores the initial prior distribution, V'^ — Vk, 
and therefore Plato has learned nothing, thus avoiding 
the paradox. Reversing the logic, in order to avoid the 
paradox, Ps must be independent of the initial state, as 
found in Eqs. ^ and 



III. DOUBLE-QUANTUM-DOT CHARGE 
QUBIT 

Consider Fig. [T] illustrating a double-quantum-dot 
(DQD) qubit, measured continuously by a symmetric 
quantum point contact (QPC). This setup has been ex- 
tensively studied in earlier papers, both theoreticallySi 
and experimentallyi^ The measurement is characterized 
by the average currents Ii and I2 corresponding to the 
qubit state |1) and |2) (the double-dot electron being in 
one dot, or the other), and by the shot noise spectral 
density Si?^ We treat the additive detector shot noise 
as a Gaussian, white, stochastic process, and assume the 
detector is in the weakly responding regime, |A/| ^ /q, 
where AJ = /i - h and /q = {h + h) /2, with QPC voU- 
age much larger than all other energy scales, so that the 
measurement process can be described by the quantum 
Bayesian formalism]^ 



A. Measurement dynamics for a non-evolving qubit 

We begin for simplicity with the assumption that there 
is no qubit Hamiltonian evolution, so that the qubit state 
evolves due to the measurement only (this can also be 
effectively done using "kicked" quantum nondemolition 
(QND) measurements^-). As was shown in Ref. [l^, at 
low temperature the QPC is an ideal quantum detector 
(which does not decohere the measured qubit), so that 
the evolution of the qubit density matrix p due to contin- 
uous measurement preserves the "murity" M while the 
diagonal matrix elements evolve according to the classi- 
cal Bayes ruleJ^ We define the electrical current through 
the QPC averaged in a time t as i{t) ~ [j^ I{t') dt']/t, 
and the quantum Bayesian equations read 



Pii(i) 

P22{t) 

M 



Pii(0)Pi(/) 



Pll(0)Pi(/)+p22(0)P2(/)' 

P22(0)P2(/) 
pil(0)Fl(/)+p22(0)F2(/)' 

/O12/VP11/O22 = const. 



(22) 

(23) 
(24) 



where the conditional (Gaussian) probability densities of 
a current / realization, given that the qubit is in |1), |2) 
are 



PiM) = ^t/^Si eM-Cl - haft/Si]. (25) 
Equations ((22|) and ((23| may be simplified by noting 

(26) 



Pllji) ^ Pll(O) 2r(t) 



P22{t) P22(0) 

where we define the dimensionless measurement result as 

r{t) = ^ m - 4] - If j\m ~ lo] dt'. (27) 

Notice that r{t) is closely related to the total charge 
passed through the QPC, and therefore r{t) accumulates 
in time. For times much longer than the "measurement 
time" Tm — '2Si/{AI)'^ (the time scale required to ob- 
tain a signal-to-noise ratio of 1), the average current / 
tends to either Ii or I2 because the probability density 
P{I) of a particular / is 



p(/)= puimii)- 



(28) 



Therefore r{t) tends to ±00, continuously collapsing the 
state to either |1) (for r 00) or |2) (for r — > —oo). 
Importantly, for the special case when the initial state is 
pure, the state remains pure during the entire process. 

This set of DQD measurement dynamics can be seen 
to be related to the general measurement formalisini^ de- 
scribed in the previous section in the following way.— For 
a fixed time t the measurement result m can be associ- 
ated with the averaged current / (or, equivalently, with 
the dimensionless quantity r). The Kraus operator Mm 
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is then diagonal in the measurement basis |1) and |2) and 
has matrix elements ^/Pi^2(X)- Notice that Pm now de- 
scribes the probability density of the result / instead of 
probability, because the measurement result becomes a 
continuous variable. 



B. Uncollapsing for the charge qubit 

In order to describe how to uncoUapse the charge qubit 
state, we note that if r{t) = at some moment t, then 
the qubit state becomes exactly the same as it was ini- 
tially, p{t) = /9(0), as follows from Eqs. (US]) and dM]). 
This of course must be the case if t = 0, i.e. before 
the measurement began, but is equally valid for some 
later time. To see why this is so from the informational 
point of view, we note that in the absence of noise, the 
measurement result from states |2) would simply be 
^i,2(i) = H/Tm- With the noise present, the measure- 
ment outcome r{t) = splits in half the difference be- 
tween states |1) and |2). Such an outcome corresponds 
to an equal statistical hkelihood of the states |1) and |2), 
and therefore provides no information about the state of 
the qubit. 

Suppose the outcome of a measurement is rg, par- 
tially collapsing the qubit state toward either state |1) 
(if To > 0), or state |2) (if < 0). The previous "no 
information" observation suggests the following strategy 
for uncollapsing: continue measuring, with the hope that 
after some time t the stochastic result of the second mea- 
surement r„(t) becomes equal to — tq, so the total result 
r{t) = tq + r„(t) is zero, and therefore the initial qubit 
state is fully restored. If this happens, the measuring 
device is immediately switched off and the uncollapsing 
procedure is successful (Fig. [T]). However r(t) may never 
cross the origin, and then the uncollapsing attempt fails. 

This strategy requires the observation of a particu- 
lar measurement result that may never materialize. The 
strategy shifts the randomness to the amount of time that 
needs to elapse in order to find the desired measurement 
result. Of course, in a given realization the measurement 
result could take on the desired value multiple times, so 
we will take as our strategy to turn off the detector the 
first time the measurement result takes on r = 0. In the 
classical stochastic physics this is known as a first pas- 
sage process, the theory of which is well developed and 
will be used below. 

In order to analyze the uncollapsing strategy perfor- 
mance, in particular to find the success probability P5, 
it is important to notice that the off-diagonal elements of 
the qubit density matrix p do not come into play when 
we consider the detector output I{t) (this is true only in 
the case of zero or QND-eliminated qubit Hamiltonian; 
we will shortly generalize to the finite qubit Hamiltonian 
case). As a result, the quantum problem can be exactly 
reduced to a classical problem by substituting pii(t) and 
P22{t) with classical probabilities, evolving in the course 
of measurement according to the classical Bayes rule. 



First "accidental" Undoing 
measurement measurement 




FIG. 1: (Color online). Illustration of the uncollapsing pro- 
cedure for the charge qubit. The slanted lines indicate the 
deterministic output of the detector in the absence of noise, if 
the qubit is in state |1) or |2). The initial measurement yields 
the result ro. The detector is again turned on, hoping that 
at some future time the measurement result r{t) = ro -\- ru{t) 
crosses the origin, at which time the detector is turned off, 
successfully erasing the information obtained in the first mea- 
surement, and restoring the initial qubit state. 



while evolution of the off-diagonal elements pi2 = P21 
can be found automatically from the murity conserva- 
tion law We therefore model the qubit by a clas- 
sical bit with probability pi = Pii(O) of being prepared 
in state "1" and probability p2 = P22(0) of being in state 
"2" . If the bit is in state "1" , the dimensionless measure- 
ment result r(t) evolves as a random walk with diffusion 
coefficient D = [Alf/^Si ^ 1/2Tm and drift veloc- 
ity vi = {AI)^/2Si = 1/Tm [see Eqs. ^ and 
For the bit state "2" the random walk of r{t) has the 
same diffusion coefficient but the opposite drift velocity 
V2 = — (A/)^/2S'/. We are given the fact that the first 
part of the measurement had the result rg (i.e. we select 
only such realizations). We need to analyze the stochas- 
tic behavior of the total measurement result r(t) during 
the second part of measurement, with most attention to 
the crossing of the zero line r{t) =0 [for convenience we 
shift t = to the beginning of the second measurement, 
so that r(0) ~ ro]. 

Let us find the probability Ps of such a crossing. We 
will first obtain it in a simple way, and then reproduce the 
result in a more complicated way, which will also allow 
us to analyze the statistics of the waiting time. For defi- 
niteness take tq > (this will be extended later). Then 
the result r{t) will necessarily cross if the bit is actually 
in the state "2", because in this case r{t = 00) — —00 
while r(t = 0) — vq > 0. The probability of being in 
the state "2" is p2 = p2e-''° / ipie''° + p2e-''°) from (j^ . 
which differs from p2 because of the Bayesian update. 
If the bit is in state "I" (this happens with probability 
Pi — 1 — P2), then r{t = 00) = -l-oo and the crossing of 
may never happen; however, it is still possible with some 
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probability Pc, which depends on rg, and also on D and 
Vi- To find Pc, let us consider an infinitesimal time step 
dt and model the diffusi on by discrete jumps in r of mag- 
nitude Ar = zL\/2D dt. After a step dt, the coordinate 
will then shift to one of two positions, r — r±, where 
r-t = ro + vidt ± \^2D dt. Each of these new coordinates 
will have its own probability of eventually crossing the 
origin, Pci'r±). Because the diffusive dynamics is gener- 
ated by choosing either r+ or r_ with equal weighting, it 
follows in the limit dt ^ that 



± ^ 



(29) 



Expanding Pc{r±) in this relation in a Taylor series, we 
find from the linear in dt term that 



DdlPc^-Vidr„Pc, 



(30) 



where 9^ and d^^ denote the first and second derivatives 
with respect to tq. Taking into account that Pc — 1 
for rp = and Pc = for ro = oo, the above dif- 
ferential equation may be easily solved to find Pc — 
exp(— uiToZ-D) — exp(— 2ro). Now collecting the prob- 
abilities of the zero line crossing for both bit states, we 
find 



Ps ^PiPc +P2 



V(pie^ 



P2e 



(31) 



The derivation for < is similar and leads to the extra 
factor e^'"" , so that the crossing probability in both cases 
can be written as Ps = e^^^°^ / {pie^° +p2e~^°). 

Thus, using the trick of reducing the quantum dynam- 
ics to the classical problem, we have found the probabil- 
ity of successful uncoUapsing for a DQD qubit with no 
Hamiltonian evolution: 



Ps 



-\ra\ 



(32) 



where pm characterizes the qubit state before the first 
measurement. We will discuss this result in the next 
subsection. Before that let us rederive it in a different 
way, using the power of the standard methods of first 
passage theoryi^ This method has recently been used to 
investigate entanglement dynamics of jointly measured 
qubit s.— 

It is convenient to scale time in units of the measure- 
ment time, r = t/TM] then the probability distributions 
(p5)) take the simple form 



A,2(r,r) ^ \/-^exp 
Zttt 



2r 



(33) 



These are the solutions of two different classical random 
walks with dimensionless drift velocity vi,2 — ±1 and 
dimensionless diffusion coefficient I? = 1/2 described by 
the Fokker-Planck equations. 



^rP^{r,T) = -V^^rP^ + DdlP^. 



(34) 



In order to solve the first passage problem, we solve first 
for the Green functions G(r, r) of the above equations 
starting from the initial condition r = tq. The solutions 
from the different drift velocities will be weighted with 
probabilities pi,2- These Green function equations are 
supplemented with an absorbing boundary condition at 
the origin (r = 0), 



G,(r = 0,t) ==0, 



(35) 



in order to account for the statistics of events that cross 
this point at least one time. Let as again start with 
assuming > and consider the other case later. The 
solution of Eq. (p4)) subject to the condition p5|) is most 
easily found by guessing: 



G,;(r,T) 



1 



exp 



exp -ViT^/b 



{r ~ ro - ViT^ 



exp 



ADt 

(r + rp - UitY 
ADt 



(36) 



In the form written above, it is obvious that the so- 
lution obeys the equation of motion and has the 
correct initial condition tq at t = (because the ab- 
sorbing boundary condition only permits r > solu- 
tions). Further inspection of the solution is facilitated 
by factoring out the free Green function, Ghcd'i',''') = 

exp[— (r — rf) — VirY / [ADt)]/ \/ AttDt to write the solu- 
tion as 



Gi{r,T) = Gfi.co(r, r) 



1 — exp 



Dt 



(37) 



One can now explicitly see that the absorbing boundary 
condition ([55]) is satisfied, and the solution is completely 
positive (as it must be to represent a probability density). 

To calculate the first passage time distribution, we first 
note that the total survival probability that the random 
walker will be in the interval r g (0, oo) at time r is given 
by -Psur(''') = /q drG{r,T). However, the only place for 
the particle to be lost from the system is at the origin. 
Therefore, the first passage time distribution P^^^ is given 
by 



£J = -drPsur = - I drdtG,{r, r). (38) 
'o 



The next step is to note that the Fokker-Planck equation 
(1341) may be rewritten as a continuity equation, drGi + 
drJi = 0. This simply means that locally, probability 
is conserved. The probability current in the continuity 
equation is Ji = —DdrGi +ViGi from p4p . Substituting 
this into ([55)1 we find the general result 



/•oo 

Fg(r) = / drdrJ, = Moo) - J.(0) - -J,(0), (39) 
Jo 

because the probability current at infinity vanishes. Ap- 
plied to our problem, we find 



P^iiir) 



V AttDt^ 



exp 



-{ro + u,T)y{ADT) . (40) 
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The probabihty Pc that the point r = is ever crossed is 
found by integrating over aU positive time to obtain 



^1 exp{-viro/D) = exp(-2ro), i ^ 1, 

I 1, i = 2. ^ ' 



This result may be understood intuitively because if the 
state is in i = 2 then the drift V2 = —1 causes r{t) to 
evolve from ro to — oo and therefore must cross at some 
time, while if the system is in state i = 1 then the drift 
?;i = 1 causes r(t) to evolve from ro to +00. Therefore, 
in order to cross r — 0, the noise term must fight against 
the drift, causing a successful crossing only occasionally. 

In order to obtain the normalized first passage distri- 
bution (conditioned on crossing), we divide (|40p by the 
probabilities (|4T|l to obtain 



ro 



: exp 



|r)V(4^r) 



(42) 

The mean first passage time may also be calculated from 
([12]) to obtain rc,i = ro/\vi\ = ro. 

Obtaining analogous results for ro < is straightfor- 
ward because the Green function for the Fokkcr-Planck 
equation (|34p is invariant under the transformation {r — > 
— r, ro — > — ro, Vi — + —Vi (or 1 ^ 2)} which is also re- 
flection symmetry about the origin. Therefore results 
(I40I41I42P can be extended using this symmetry. Com- 
bining results, we can now calculate the total uncoUaps- 
ing probability, Ps = PiPc,i +P2Pc,2 to obtain the result 
P2p in this new, more powerful way. 



C. Discussion 

We now discuss the physical meaning of the result 
When the first measurement result indicates a particular 
qubit state with good confidence (|ro| ^ 1), the probabil- 
ity of success Ps given by Eq. (|32p becomes very small, 
eventually becoming P5 = for a projective measure- 
ment, realized for ro = ±00. This recovers the traditional 
statement of irreversibility. In the other limit of rg = 0, 
the success probability is unity because no time needs to 
elapse - the state is already undisturbed. We stress that 
the possibility of uncoUapsing as well as our formalism 
requires a quantum-limited detector, i.e. one that intro- 
duces no additional dephasing to the system. For such a 
detector measuring a pure state, the state remains pure 
throughout the partial collapse, and the uncollapse. We 
also note that if the qubit is entangled with other qubits, 
the uncoUapsing restores the state of the whole system. 

Let us compare the general upper bound (|13p for the 
success probability Ps with the result ([5^ . Substitut- 
ing the probabilities in the bound with probability 
densities, we find 



Ps < 



min{Pi(/),P2(/)} 

Pl{I)Pin,ll + P2{I)Pin,22 '' 



(43) 



where / corresponds to the measurement result rp. The 
bound coincides with Eq. ([5^ because Pi{I)/P2{I) — 
g2ro^ This means that the "wait and stop" strategy an- 
alyzed above is optimal in the sense that it reaches the 
upper bound. 

It is also instructive to not specify the result of the first 
measurement, but to find the total probability p^°^^^ [see 
Eq. (PO)) ] that the initial qubit state can be restored after 
a measurement for time t. This is given by averaging 
Ps in Eq. ([5^ over the results ro with the corresponding 
weights (1^51) . This averaging is technically easier using 
the form of Eq. and gives 



P. 



-total 



d/min{Pi(/),P2(/)} = 1 - erf 



2T, 



M 



(44) 

which depends only on the "strength" t/Tm of the first 
measurement, but not on the initial state, as expected 
from the discussion in Sec. [TTdI Notice that the resuh 
(|44|) reaches the upper bound ((20|l because ([32|l reaches 
the upper bound ((Tc 



D. Statistical features of the time to uncollapse 

In addition to the probability of success, the complete 
solution of the first-passage problem given above now al- 
lows us to specify further information about the uncol- 
lapsing process. In particular, an important question for 
an experimental implementation of this idea is how long 
it is necessary to wait. 

Since the distribution of the first passage time ((42|) 
does not depend on the bit state, it directly gives the 
distribution of the waiting time to uncollapse any qubit 
state. Therefore, rescaling back the time axis in Eq. (I42p , 
we find the waiting time distribution is 



tit) 



exp 



<\ro\^t/TMr 



2t/T, 



M 



(45) 



This distribution is normalized, since we consider only 
successful attempts of uncoUapsing. The fact that the 
distribution is independent of the initial qubit state is 
not surprising, since otherwise a successful uncoUapsing 
instance would give us an information about the qubit 
state (see discussion in Sec. IIIPp . 

Using the distribution (j45p , we can find the mean wait- 
ing time to uncollapse 



Twait = Tm \ro\ 
the standard deviation 



(46) 



ATwait =Tm/N, (47) 
and the most likely waiting time (which maximizes Pwait) 



M 



9/4 -3/2 



(48) 



10 




Time (T|^) 

FIG. 2: (Color online). Probability distribution of the time 
required to undo the measurement. Different plots are for 
different values of ro. 

The distribution of the waiting times is plotted in 
Fig. [2] for several values of tq. Note that it has a long 
tail, which makes the average value Twait to be longer 
than the most likely value T/. 



E. Evolving charge qubit 

Let us now extend the first example of a single charge 
qubit measured by QPC, by including internal evolution 
of the qubit via a qubit Hamiltonian, 

HQB = -{el2)a, + Ha.,, (49) 

where e is the energy asymmetry between the quantum 
dot levels, and H is the tunnel coupling between the dots. 
In this case Eqs. ([^^ - ([M]) are no longer valid and should 
be replaced by the Bayesian equations^* (in Stratonovich 
forniS^) 

2A/ 

Pii = -P22 = -2Hlm pi2 + P11P22 [I{t) - 4] ,(50) 

01 

P12 = iepi2 + iH[pii - P22) 

-{Pll- P22)^[I{t)-Io\pi2, (51) 

where I{t) is the QPC current, 

i{t)^pii{t)h+P22{t)h+m. (52) 

containing white noise S^{t) with spectral density Si, and 
we use h — 1. These evolution equations are nonlinear 
and not very simple to deal with. To discuss the undoing 
of a continuous measurement, it is more convenient to 
use a non-normalized density matrix a, which has an 
advantage of dealing with linear equations. 



We rewrite Eqs. ((50)) - (|5T1) in the form 

p = cr/Trcr, (53) 
an = -2HliTvar2 - (Tii^[I{t) - hf , (54) 

cT22 = 2Hlmai2-<J22^[I{t)-l2f, (55) 
(T\2 — iea\2 + ill(a\\ — 022) 

so that cr(0) = p{0), while the ratio a{t)/p{t) decreases 
with time and is equal to the normalized probability den- 
sity of the corresponding realization of the detector out- 
put I{t'), < t' < tM In the language of general quan- 
tum measurement this formulation corresponds to omit- 
ting the denominator in Eq. ([9]). Notice that we still con- 
sider an ideal detector, so an initially pure state remains 
pure, |cri2p = aii(T22-^° 

A casual inspection of Eqs. ([54l) - ((56l) shows that they 
are seemingly not well-defined because the terms [I{t) ~ 
-fo.1,2]^ contain the term ^(t)^ — 00 [from the relation 
(C(*')C(0)) = {Si/2)6{t)]. This divergence is artificial be- 
cause there will always be a small correlation time T 
of the noise and/or a finite detector bandwidth B (corre- 
sponding to T = 1/45), so there will be a large but finite 
constant C = {^{t)"^) — Sj/AT contained in terms of the 
form [/(t)-/o,i,2]^- It is easy to see that Eqs. (I54l)-(l56|) do 
not change if we subtract the same constant from these 
terms [I{t) - Iox2? ^ [I{t) - Io,i,2? - C. This can be 
shown by considering another unnormalized density ma- 
trix rj = crexp(t/r). Writing the linear Bayesian equa- 
tions l[54 |l -([56 l) in the form &ij — /y [a-], the equations 
transform to 7/y = [77 exp(— i/T)] exp(t/T) -I- ryy/T un- 
der the change of variables. The unnormalized Bayesian 
equations are linear in the density matrix elements (Jy , 
so the exponential factors cancel out. The new equa- 
tions are thus the same as the old ones with a constant 
C = S'//4T subtracted from the [/(t)-/o,i,2]^ terms. The 
unspecified constant T in the density matrix transforma- 
tion may be chosen to be the short correlation time T 
discussed above, thus canceling the large term and mak- 
ing Eqs. (|54)) - (|56l) well-defined. The only price to be paid 
for this transformation is an altered normalization, that 
will cancel in the normalized density matrix (|53p . 

For a particular realization of the detector output 
lit'), < t' < t, Eqs. dMl-dSSl) define a linear map 
cr(0) o'(i), corresponding to a particular Kraus op- 
erator Mm (which, therefore, can be denoted as M^jy). 
For the uncoUapsing we have to realize the map, corre- 
sponding to the inverse Kraus operator CM^| (see Sec. 
ITT|) . It is obvious that in contrast to the case of the non- 
evolving qubit, this cannot be done by simply continuing 
the measurement and waiting for a specific result. The 
reason is that now the map is characterized by 6 real 
parameters (8 parameters for a linear operator CM|^j 
with neglected overall phase and normalization) , instead 
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of 1 parameter for the non-evolving case [see Eq. p4)) and 
(PS)) ] . We will discuss a little later how the 6-parameter 
uncoUapsing procedure can be realized explicitly. Before 
that we discuss how to find the operator Mrn in a more 
straightforward way, than from Eqs. (j54[) - (|56p . 

Let us consider only the evolution of (unnormalized) 
pure states \il;{t)) = a{t)\\) ^ I5{t)\2) , so that cr = 
Then Eqs. can be rewritten as 

a = +j£a-iF/J-«^[/(t)-/i]2, (57) 



/3 = -z|/3-ii/a-/3^[/(t) 



(58) 



where the infinite part of I'^it) can be canceled in the 
same way as discussed above. The linearity of these equa- 
tions guarantees that for any given realization of lit), 
it is sufficient to solve (|57l58p for the initial states |1) 
and 1 2) in order to find the solution for an arbitrary ini- 
tial state of the qubit. Defining v\ — a\{t)\\) + /3i(t)|2) 
as the solution of (|57l58p for initial state |1), and V2 = 
a2(t)\\) /32(i)|2) as the solution of (|57l58p for initial 
state 1 2), we can write the solution for an arbitrary initial 
state IV'm) = = a\\) -I- 6|2) as \i;{t)) = avx + hv-i- 

Therefore, the Kraus operator M^/j for a given realiza- 
tion of /(t), in the |1), |2) basis is 



Q!l(t) 



a2{t) 



(59) 



For the uncoUapsing we need to apply the Kraus op- 
erator CM^j^y, which maps the state vi onto C|l) and 
the state V2 onto C|2). The reason why we need a non- 
unitary transformation is that the vectors vi and V2 are 
in general non-orthogonal and have different norms. Ge- 
ometrically, such a transformation can be done by us- 
ing a relative shrinking or stretching of two orthogonal 
axes (found for a given M^jy), which would make vi and 
V2 orthogonal and equal in norm, followed by a unitary 
transformation (this would correspond to the decompo- 
sition of the form U^fE - see Sec. Ill B|) . However, for 
a practical realization of uncoUapsing it is most natu- 
ral to use the shrinking or stretching of the axes |1) and 
1 2) via a continuous QND measurement with the QPC 
in the way considered above for a non-evolving qubit. 
In this case the uncoUapsing procedure can be done in 
three steps (see Sec. Ill B| ): unitary evolution V , continu- 
ous QND measurement (where the qubit Hamiltonian is 
turned off, e = H = 0) described by a diagonal matrix 
i, and a final unitary operation U . (In the notation of 
Sec. Ill B| V corresponds to V}^U^, and U corresponds to 
U] .) These operators should satisfy 



ULV = CM^^/p 



(60) 



and it is easy to find U , L, and V explicitly by recognizing 
Eq. (|60|) as a singular value decomposition of the operator 
CAf^j (recall here that L is diagonal; also notice that 



the standard form for the singular value decomposition 
is slightly different, with V denoted as V^). 
To find L explicitly, we notice that 



(61) 



which is simply the diagonalization of C ^Af|^jAf{/].. 
Therefore, 



L = C 



where 



\z'^' 



a; 



-1/2 



L = cr'+ 



x-'^' 



A 



-1/2 



(62) 




are the eigenvalues of the operator M^j^M^jy and the 
vectors Vi are defined above Eq. ([55]) . The Cauchy- 
Schwartz inequality, \vi ■ v^\'^ < \\vi\\'^\\v2\\^ , guarantees 
the non-negativity of A_. (The notation vi ■ is used 
for the inner product {v2\vi) of V2 and vi). 

To find U, we use Eq. again and see that 

the columns of U are composed of the eigenvectors of 
M^jyM^jj [the sequence of columns depends on the 
choice in Eq. (l62|) ]. Finally, V is given hy V = 
WCL-'^M^/y For brevity we will not show the matrices 
U and V explicitly. 

In the physical realization of the uncoUapsing proce- 
dure the measurement step L can be performed in ex- 
actly the same way as in Sec. IIII Bl Comparing Eq. ((62|) 
with Eqs. (|26p and ([M]) . we see that the continuous mea- 
surement by the QPC should be stopped when the di- 
mensionless measurement result r(t) reaches the value 



ri = In v/A+/A_ > or ra = In \/X-/X+ < 0, (64) 

for the first and second choice in Eq. (p^ . respectively 
(the choice should be made beforehand, since it deter- 
mines operation V). As previously mentioned, the con- 
stant C is not important here because the physical state 
is always normalized. The procedure fails if the desired 
result is not reached during the continuous measurement. 

The unitary operations V and U can be practically re- 
alized in three substeps each: z-rotation on Bloch sphere 
by applying non-zero energy asymmetry e for some time, 
y-rotation by applying non-zero tunneling H, and then 
one more z-rotation. However, the last z-rotation of V 
and the first z-rotation of U are simply added to each 
other (since L does not change the relative phase of the 
state components or, equivalently, the azimuth angle on 
the Bloch sphere). The corresponding trivial degree of 
freedom can be eliminated, for example, by realizing the 
operation V in only two substeps, without the second 
z-rotation. 

Let us count the number of real parameters, character- 
izing the uncoUapsing procedure. Since V and U together 
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provide 2x3 — 1 = 5 parameters, and the desired result 
r in the measurement step adds one more parameter, 
the overall number of parameters is 6. As expected, this 
is exactly the needed number of parameters characteriz- 
ing an arbitrary Kraus operator for the qubit (neglecting 
normalization and overall phase). Let us also mention 
the fact from linear algebra that the singular value de- 
composition (j60p is unique in the non-degenerate case 
(A+ > A_ > 0), up to the permutation of singular val- 
ues [corresponding to the choice in Eq. ([S^ ] and arbi- 
trary phase factors in columns of U, with compensating 
changes in V (this corresponds to the discussed above 
compensation of z-rotations). 

Now let us discuss the probability Ps of the successful 
uncoUapsing. From the general theory discussed in sec- 
tion II [Eq. ([5])], it is bounded from above by a fraction, 
Ps ^ ^^'^ P{i} / P{i}{\'fpin)) , in which the denominator is 
the probability density of the given realization I{t) for 
the initial state \^p,n) = a\l)+b\2) (with \a\'^ + \b\'^ = 1), 
while the numerator is this probability minimized over 
all initial states. So, the denominator is given by the 
squared norm of the final state \tp{t)) = avi + hv2, 



P{I}{\i^^n)) = \\avi+bv2\?. 



(65) 



while the numerator is given by minimizing (j65p over 
all normalized initial states. It is easy to see that this 
minimum is equal to the minimum eigenvalue A_ of the 
operator m|^jM{/}, given by Eq. ((63| : therefore. 



Ps < 



A- 



\\avi + bv2\\^' 



(66) 



Converting this result into the language of density ma- 
trices and simultaneously generalizing it to an arbitrary 
initial state pin, we obtain the bound 



Ps < 



/3m,ll||wi||2 -I- pm,22||w2||^ + 2Re[p^ 



'm,12 Vi ■V2\ 



(67) 



in which the numerator is the explicit expression 
for A_. It is easy to check that this result reduces to 
the bound (I43p in the non-evolving case, in which vi = 

/fUO,0)^ and V2 = (O, y^lWff . 

The uncoUapsing procedure discussed in this subsec- 
tion is optimal in the sense that it corresponds to the 
upper bound of Eq. ([S7|) . To prove this statement, in- 
stead of calculating Ps explicitly, let us use the fact (see 
Sec. Ill Pp that the product PsP{i} cannot depend on the 
initial state. Therefore, it is sufficient to prove the op- 
timality of Ps only for one initial state. Let us choose 
the state [■0i„) that is the eigenvector of M'^^j^M^^iy, cor- 
responding to the eigenvalue A_. Then after the first 
measurement (operator M^/j.) and the unitary operation 
V it is transformed into one of the basis states [|1) or |2) 
for the first or second choice in (|62p. respectively]. Re- 
call that for the QND (non-evolving) measurement case. 



r{t) — >■ 00 for the initial state |1), while r{t) —oo for 
the initial state |2). The crossing thresholds ([M]) indicate 
that the measurement L is always successful because r{t) 
necessarily crosses the desired value (which is positive for 
|1) and negative for |2), as discussed above). Therefore, 
the uncoUapsing success probability for this special state 
is 100%, that is equal to the upper bound ([S7)l . As men- 
tioned above, the optiniality of the procedure for this 
special state also proves the optimality for any initial 
state. 



Obviously, an uncoUapsing procedure can also be 
non-optimal. As an example, let us consider a pro- 
cedure which realizes the desired mapping {wi,W2} — >■ 
{C|l), C|2)} using two measurements instead of one. The 
goal of the first measurement is to map {vi,V2} into an 
orthogonal pair of vectors, while the goal of the second 
measurement is to equalize their norms, keeping them 
orthogonal. The first goal can be achieved by stretch- 
ing/shrinking of the Hubert space along any axis u of the 
form vi+cv2{vi-V2)/\vi''V2 \ with an arbitrary positive real 
number c (it is easy to visualize this procedure of mak- 
ing two vectors orthogonal by assuming the space of real 
vectors, for which the axis u is geometrically in between 
vi and V2'-, the same geometrical idea works for complex 
vectors). We recall that measurement for a non-evolving 
qubit (Sec. |TTTX| stretches (squeezes) the |1) axis, while 
squeezing (stretching) the |2) axis. Therefore, the first 
goal can be achieved by a unitary operation which ro- 
tates u into |1), followed by a continuous measurement 
(with a QPC) of a non-evolving qubit, to be stopped 
when the mapped vectors become orthogonal. After the 
vectors {vi,V2} are transformed into an orthogonal pair 
by the (successful) first measurement, the second part of 
the procedure should stretch/shrink the 2D Hilbert space 
along the resulting vectors to make them equal in norm. 
This can be done similarly, by a unitary rotation and 
partial measurement of a non-evolving qubit. Finally, 
another unitary operation can be used to map the re- 
sulting pair of vectors into {C|l), C|2)}, thus completing 
the uncoUapsing procedure. Notice that both measure- 
ments are performed in the "wait and stop" manner, and 
both measurements should be successful to realize the 
uncoUapsing. While the successfully uncoUapsed state is 
stUl perfect in this procedure, the probability of success is 
lower than the bound (|67p . To prove this non-optimality, 
let us again use the initial eigenstate [V'm): which corre- 
sponds to eigenvalue A_, so that the bound (|67p is 100%. 
Then the probability of success for the first measurement 
is in general less than 100% (it is 100% only for one spe- 
cific axis discussed previously, while here we consider a 
range of possible axes by allowing c to vary). Thus, the 
success probability is less than 100% for this special state, 
and therefore Ps is below the bound (j67p for any initial 
state. 
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FIG. 3: (Color online), (a) Schematic of a phase qubit con- 
trolled by an external flixx (j>e and inductively coupled to the 
detector SQUID, (b) Energy profile V{(j)) with quantized lev- 
els representing qubit states. The tunneling event through 
the barrier is sensed by the SQUID. 



IV. PHASE QUBIT 

The second explicit example of erasing information and 
uncoUapsing the wavefunction is for a superconducting 
phase qubit. The system is comprised of a supercon- 
ducting loop interrupted by one Josephson junction [Fig. 
[3]Ja)], which is controUed by an external flux (jje in the 
loop. Two qubit states |1) and |2) [Fig. [3^b)] correspond 
to two lowest states in the quantum well for the poten- 
tial energy V{(j)) where (j) is the superconducting phase 
difference across the junction. (Notice that the standard 
notation for the phase qubit states is |0) and |1);'^^ how- 
ever, we use |1) and |2) for consistency with the previous 
section.) The qubit is measured by lowering the barrier 
(which depends on 0e), so that the upper state |2) tun- 
nels into the continuum with the rate F, while the state 
|1) does not tunnel out. The tunneling event is sensed 
by a two-junction detector SQUID inductively coupled 
to the qubit [Fig. [3l^a)]. Transitions between the levels 
|1) and 1 2) can be induced by applying microwave pulses 
that are resonant with the energy level difference. 



A. Partial collapse 

For sufficiently long tunneling time t, Tt ^ 1, 
the measurement is a (partially destructive) projective 
measurement!^ the system is destroyed if the tunneling 
occurs, while if there is no record of tunneling, then the 
state is projected onto the lower state |1). This mea- 
surement technique is remarkable in the fact that the 
wavefunction is collapsed if nothing happens. A more 
subtle situation arises if the barrier is raised after a finite 
time t ^ F^^: then the measurement is only partial and 
therefore is of a POVM-type>^ The system is still de- 
stroyed if tunneling happens, while in the case of no tun- 
neling (which we refer to as a null-result measurement) 
the state is partially collapsed. This situation may be 
described with a two-outcome POVM, with elements £"„ 
and Ey, where n denotes the null result, and y denotes 
the affirmative (tunneling) result. The POVM elements. 



given in the |1), |2) basis are 




with the obvious completeness relation £"„ + Ey ^ 1. 

It is interesting to notice that while En = Mj^Mn cor- 
responds to a Kraus operator M„ (see discussion below) , 
no meaningful Kraus operator My can be introduced for 
the POVM element Ey, because in the case of a tunnel- 
ing event the system leaves its two-dimensional Hilbert 
space and becomes incoherent (so that a single Kraus op- 
erator cannot be introduced even in the extended Hilbert 
space). However, this is not important for us because we 
are interested in the null-result case only. 

Limiting the unitary operation in the decomposition 
^ to be the phase factor only, we may expect the 
null-result Kraus operator Af„ to be of the form Af„ — 
diag{l, e~'"*/^e^*'^}. In a simple modeP**'?- the phase 
(p is zero in the rotating frame, which compensates for 
the energy difference of the states |1) and |2). In a real 
experiment,"^? however, this energy difference changes in 
the process of measurement because it is affected by (j)e, 
and therefore even in the rotating frame the phase tp 
is non-zero. Correspondingly, the qubit density matrix 
changes after the null-result measurement as 

Plljt) _ PinM Pl2{t) _ e^'^^'Vro,12 

P22(<) e-rtp-„ 22' y/pn{t)p22{t) ^yPin,llPin,22 ' 

(69) 

In the real experiment'^^ the situation is even more com- 
plex because the tunneling rate gradually changes in 
time; also, instead of controlling the measurement time 
t, it is much easier to control the tunneling rate. As a 
result, the measurement should be characterized by the 
overall strength = l — exp[— J*T{t')dt']. Nevertheless, 
for simplicity, we use here the physically transparent lan- 
guage of Eq. with e~^* understood as 1 — pt. 

Up to such changes of notation, the coherent non- 
unitary evolution (j69p has been experimentally verified 
in Ref. [l^ using tomography of the post-measurement 
state. The state tomography consisted of 3 types of ro- 
tations of the qubit Bloch sphere, followed by complete 
(projective) measurement. Actually, in the experiment 
it was not possible to select only the null-result cases, 
because it was not possible to distinguish if a tunneling 
event happened during measurement or during tomogra- 
phy. However, a simple trick of comparing the protocols 
with and without tomography made it possible to sepa- 
rate the null-result cases. 

Notice that except for the effect of extra phase (p{t), 
the qubit evolution (j69p is similar to the qubit evolution 
in the example of Sec. IIII Al in particular, it also rep- 
resents an ideal measurement which does not decohere 
the qubit. Formally, the evolution (|69p corresponds to 
the measurement result r = Ft/2 in Eq. As will be 

shown later, the probability to undo the measurement is 
still given by Eq. using this value of r. 
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B. Uncollapsing 

We will now describe how to undo the state distur- 
bance ([nH) caused by the partial collapse resulting from 
the null-result measurement. The undoing of this mea- 
surement consists of three steps: (i) Exchange the am- 
plitudes for the states |1) and |2) by application of a mi- 
crowave TT-pulse, (a) perform another measurement by 
lowering the barrier, identical to the first measurement, 
(Hi) apply a second 7r-pulse. If the tunneling did not hap- 
pen during the second measurement, then the informa- 
tion about the initial qubit state is canceled (both basis 
states have equal likelihood for two null-result measure- 
ments). Correspondingly, according to Eq. (j69p (which 
is applied for the second time with exchanged indices 
1 <-> 2), any initial qubit state is fully restored. An added 
benefit to this strategy is that the phase ip is also can- 
celed automatically; the physics of this phase cancellation 
is the same as in the spin-echo technique for qubits. 

It is easy to mistake the above pulse-sequence as simply 
the well known spin-echo technique alone. We stress that 
this is not the case: Spin-echo deterministically reverses 
an unknown unitary transformation (arising usually from 
a slowly varying magnetic field) without gaining or losing 
any information about what that state is. Our strategy 
is probabilistic and requires erasing the classical infor- 
mation that one extracts from the system to begin with. 
It is a (probabilistic) reversal of a known non-unitary 
transformation - and therefore quite different from spin 
echo. 

The success probability Ps for the uncollapsing strat- 
egy described above may be calculated by noting that 
it is just the probability that the second measurement 
gives a null result. If we start with the qubit state pi„, 
the state after the first measurement is given by Eq. 
and after the 7r-pulse the occupation of the upper state 
is P22 = Pin, 11 /[pin, 11 + Pin, 22 e^'"*]. The success proba- 
bility is simply the probability that the second tunneling 
will not occur, Ps = pn + P22e^^* = 1 - ^22(1 - e"'"*), 
which can be expressed as 

Ps = , (70) 

Pm,ll + e Pm, 22 

and formally coincides with Eq. ([5^ for r = rt/2. We 
can verify that this strategy is optimal by using the £"„ 
POVM clement jMl) together with the general result (fTa]) 
to find the upper bound for the success probability P5. 
The numerator of is the smallest eigenvalue of En, 
which is exp(— rt), while the denominator is Tri?„pi„ = 
Pm,ii + e~^^pin.22, giving a Ps that coincides with Eq. 
(|70p . and thus confirming the optimality of the analyzed 
uncollapsing procedure. 

The total uncollapsing probability Ps of two null re- 
sults (see Sec. II D) is 

Ps = {p^nM + e-^' p^n,22)Ps = e"^*. (71) 

As expected (see Sec. II D) this probability does not de- 
pend on the initial state. 



Uncollapsing of the phase qubit state has recently been 
experimentally realized by Nadav Katz and colleagues in 
the lab of John Martinis, at UC Santa Barbaraf^i The 
experimental protocol was slightly shorter than that de- 
scribed above: it was missing the second 7r-pulsc, so the 
uncoUapsed state was actually the 7r-rotation of the ini- 
tial state. Shortening of the protocol helped in decreas- 
ing the duration of the pulse sequence, which was about 
45 ns, including the state tomography. Since the qubit 
energy relaxation and dephasing times were significantly 
longer, Ti — 450 ns and Tj* = 350 ns, the simple the- 
ory described above was sufficiently accurate. The same 
trick as for the partial-collapse experiment^ was used to 
separate tunneling events during the first, second, and 
tomography measurements, because the detector SQUID 
was too slow to distinguish them directly. 



The uncoUapse procedure should restore any initial 
state. However, instead of examining all initial states 
to check this fact, it is sufficient to choose 4 initial states 
with linearly independent density matrices and use the 
linearity of quantum operationsj^^ In the experimentii 
the uncollapse procedure was applied to the initial states 
(|1> + |2))/\/2, (|1) - i|2))/\/2, |1), and |2), and then 
the results were expressed via the language of the quan- 
tum process tomographyi^i^S, (QPT). The experimentalii 
QPT fidelity of the uncollapsing procedure was above 
70% for pt < 0.6. A significant decrease of the uncol- 
lapsing fidelity for larger measurement strength pt, es- 
pecially for Pt > 0.8, was due to finite Ti time and the 
fact that the null-result selection preferentially selects the 
cases with energy relaxation events, so that the proce- 
dure should no longer work well when 1 — pt becomes 
comparable to the probability of energy relaxation. 



As mentioned above, the uncollapsing procedure de- 
scribed in this subsection is theoretically optimal in the 
sense that it maximizes the bound (|13p for the uncol- 
lapsing probability. An example of a non-optimal uncol- 
lapsing for a phase qubit was considered in Ref. [H. It 
was shown that if the measurement process is performed 
simultaneously with Rabi oscillations, then in the null- 
result case the initial state is periodically restored. The 
non-optimality of uncollapsing for such a procedure is due 
to measurement of an evolving qubit, which corresponds 
to a sequence of many measurements; a similar reason 
for the non-optimality of the two-step uncollapsing was 
discussed at the end of Sec. IIII El 



Exact uncollapsing requires an ideal detection, which 
does not decohere a quantum state; Eq. ([Mjl corresponds 
to such an ideal detection. However, if various decoher- 
ence mechanisms are taken into account;^ then only im- 
perfect uncollapsing is possible. The theory of imperfect 
uncollapsing is a subject of further research. 
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V. GENERAL PROCEDURE FOR ENTANGLED 
CHARGE QUBITS 

Let us present an explicit procedure which can be used 
in principle to undo an arbitrary measurement M„ of 
any number N of entangled qubits with maximum prob- 
ability. For simplicity we consider double-quantum-dot 
charge qubits and assume that any unitary transforma- 
tion can be used in the procedure. If the operator Mm 
was produced by a one-qubit measurement, and other 
entangled qubits were not experiencing a Hamiltonian 
evolution, then the formalism of Sec. IIII Al is essentially 
unchanged,— and uncoUapsing of the measured qubit 
leads to the restoration of the whole entangled state. If 
the operator Af,„ was produced by a one-qubit measure- 
ment, while other qubits were evolving in a unitary way 
but not interacting with the measured qubit, then uncol- 
lapsing is also easy: we should uncoUapse the measured 
qubit in the usual way (Sec. Ill) and should apply inverse 
unitary transformation for other qubits. In this section, 
however, we do not consider these simple special cases; 
the goal is to undo an arbitrary Kraus operator Mm- 

Let us decompose Mm as Aim = UmV^m [see Eq. 
([3])]. Reversing the unitary operation Um can be done 
in the regular Hamiltonian way, so the nontrivial part 
is undoing the operator. We recall the diago- 

nalization of Em is given by Em = with 
vectors \i) forming an orthonormal basis. As discussed 
in Sec. II B, for the optimal uncoUapsing which maxi- 
mizes the success probability, we have to perform a pro- 
cedure corresponding to another measurement operator 

-1/2 



L — y miuj p^™-* Em'''^ which is also diagonal in the basis 

\i) with corresponding matrix elements La — = 

{jaiUj p^p^ ) /pf^'^ , all of which are between and 1. 

Given N qubits, i ranges from 1 to 2^. Notice that L is 
obviously Hermitian. 

Our procedure is to realize L with a sequence of null- 
result measurements and unitary operations. Shown in 
Fig. [4] is an illustration of the physical set-up that is used 
for the measurements: a QPC (tunnel junction) capac- 
itively coupled to non-evolving DQD charge qubits. 
We assume that the QPC is tuned to a highly nonlin- 
ear regime, for which no electron can tunnel across the 
QPC barrier on experimentally relevant time-scales un- 
less all qubits are in the state |1). We name this multi- 
qubit state |1) = |1, 1, . . . , 1). Such a regime is possible 
because of the exponential dependence of the tunneling 
rate on QPC barrier height, while the barrier height de- 
pends linearly on the states of the coupled qubits. Of 
course, this regime is not quite realistic; however, we dis- 
cuss the procedure in principle. We also assume that 
even for the A^-qubit state the rate 7 of electron 
tunneling through the QPC is rather low, so that we 
can distinguish single tunneling events (technically, this 
would require an additional single-electron transistor). 
If we perform the measurement during time t and see no 
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FIG. 4: (Color online.) Schematic set-up for uncoUapsing of 
A*' entangled qubits. The tunnel junction detector (QPC) is 
in a strongly nonlinear regime, so that an electron can tunnel 
through it with rate 7 only when all qubits are in state 



tunneling through the QPC, then similarly to the case of 
Sec. IV A, the corresponding null- result Kraus operator 
M„ shrinks the |1) axis of the Hilbert space by the factor 
e 1'*/^, while leaving all perpendicular axes unchanged. 
For the matrix elements this means (1|M„|1) = e"'''*/^, 
(V^^|M„|^jL^ = (V,^|M„|1) = (l|Af„|^/) - 0, 

where we introduced a set of 2^ — 1 states |?/;^) spanning 
the subspace orthogonal to 

The general strategy to implement the operator L is 
the following. We first note that in the basis \i) that 
diagonalizes L, this diagonal matrix can be represented 
as a product of 2^ diagonal matrices, where each term 
in the product has all diagonal entries as 1, except the 
ith entry: diag{l, 1, . . . , Lu^ . . . , 1}. Each of these ma- 
trices may be interpreted as a separate Kraus operator 
that can be sequentially implemented. Thus, the explicit 
physical procedure consists of 2^ steps, each of which has 
3 substeps. First, we apply a unitary transformation Ui 
which transforms the first basis vector |i = 1) into the 
state Then the evolution of all qubits is stopped, 
and the detector is turned on for a time ti. This time 
is chosen so that the null-result Kraus operator Lf^^^ has 
the desired matrix element = Ln; this condi- 

tion yields ti = — 27^^1nLii. The measurement is then 
followed by the reverse unitary, u\,to take the state |1) 
back to state \i = 1). This 3-substep procedure is then 
repeated for i — 2, 3, ... , 2^, sequentially transforming 
the state \i) to |1) with unitary C/^, and performing mea- 
surement with the detector for a time ti = —2j~^ In La, 
followed by the reverse unitary, uj . This sequence of 
steps decomposes the uncoUapsing operator L as 



(72) 



The uncoUapsing procedure is successful only if there 
were no tunneling events in the QPC. By construction, 
the success probability P5 for this procedure maximizes 
the general bound (fT5|) . 

The success probability P5 for the uncoUapsing process 



with Pin = LpLyTT{L^Lp) and p = Uj^pmUm, 
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can be calculated as 

Ps = Tr(LtZp) ^ ^ Zip,, = ^ p,, exp(-7i,), (73) 

i i 

where pij are the matrix elements of p in the basis 
which diagonalizes L. We may also find this result from 
another perspective by realizing that the success proba- 
bility is simply the product of the null-result probabilities 
p'g of all 2^ measurements, 

2" ~ -7t- I Y^2" 

P - TT r,(*) - ^j=^P"^ +l^J=t+lP0i 

- liPs , Ps - i-1 , — — — ' 

(74) 

where the expression for p^g comes from comparing the 
traces of unnormalized density matrices after each of 2^ 
steps of the procedure. It is instructive to show explic- 

itly that this expression for is equal to the expected 
expression 

p« = l-(l-e-^*0/5if, (75) 

in which p^*^ is the normalized density matrix before the 
ith step of the procedure (after i — 1 null-result steps). 
This can be done if we prove the relation 

k k 

- (^{l - e-^*-)] = 1 - E(l - e-^*0/5.., (76) 

i=l i=l 

[notice that the right-hand-side of this equation is equal 
to the numerator in Eq. (ffi]) with substitution k ^ i]. 
Equation ((76)) can be proven by induction using the re- 
lation pf^ = Pii/n}=i[l - (1 - e~''*OPjj-^]' which can be 
easily derived recursively, pf^ Pii^^\ starting from 
p"^^ — Pa. In this way we show consistency between the 
null-result probabilities given by Eqs. ([Ti]) and ([75]) . per- 
mitting the calculation of Ps in two independent ways. 

Let us mention again that the uncoUapsing procedure 
considered in this section reaches the upper bound 
for the success probability Pg, that can be seen both by 
construction and explicitly. 

VI. RECENT DEVELOPMENTS IN 
WAVEFUNCTION UNCOLLAPSE 

Before concluding, we wish to give a summary of some 
interesting recent developments in this area of research. 
We will briefly discuss one theory proposal and one ex- 
periment. 

A. Spin qubit 

The examples given above mainly concern quantum 
dot charge qubits. It is a natural question if a similar 



kind of partial coUapse/uncoUapse can be carried over to 
spin qubits. An analysis of this situation was carried out 
by Trauzettel, Burkard, and one of the authorsi^ There 
it was shown how an uncollapse measurement can be re- 
alized using a scheme similar to the recent experiments 
by Koppens et al.'^^ The essential idea of the spin-qubit 
experiments^i^i^^ji^ is to manipulate and measure the 
spin of a single electron through the charge degree of free- 
dom. This technique circumvents the otherwise difficult 
problem of controlling the weakly interacting spin. While 
we refer the reader to Refs. 39,40,41,42,43 for the details, 
we will give a simplified thumb-nail sketch of the physics 
here. 

The qubit is encoded with two electron spins, where 
each electron is confined in a separate quantum dot. In 
contrast to our charge qubit discussion, these dots are 
open, with electrons able to enter and leave. Electrical 
bias is applied across this double quantum dot leading 
to charge transport. Electrons can tunnel sequentially, 
but spin blockade restricts transport to situations where 
the two electrons form a spin singlet (0, 2)5' on the right 
dot while the spin triplet (0, 2)T is outside the transport 
energy window due to the large single quantum dot ex- 
change energy [here (n, m) refers to n electrons on the 
left dot and m electrons on the right dot] . This blockade 
physics provides an interesting initialization procedure of 
the quantum register - when the single-electron current 
stops fiowing, we are confident that the two-electron state 
is in a (1, 1)T state, because in the absence of spin flip 
processes, the tunneling transition to the (0, 2) state is 
forbidden. From this conflguration, it is possible to ma- 
nipulate the system by applying electron spin resonance 
pulses,^-2 transitioning the state to have overlap with the 
singlet state. Thus, the electron on the left dot may tun- 
nel (with rate F) to the right dot and exit the system, 
giving rise to a small electrical current at the drain when 
this process is repeated many times. Of course, this will 
happen with some probability controlled by the overlap 
of the state with the singlet. 

Drawing on our experience with the phase qubit, it 
is clear how to devise a weak measurement experiment 
and an uncoUapsing experiment: the allowed transition 
can be permitted for a time of one's choosing and then 
forbidden by detuning the energy levels with a voltage 
pulse to one of the quantum dot's gates. In this way 
one can weakly probe the two-electron state, and in the 
null-result case (no single electron tunneling) partially 
collapse it to the triplet subspace. In order to propose 
the uncollapse part of the experiment, it is easiest to 
consider the case when the nuclear spins quickly ad- 
mixed the singlet state with a triplet state, permitting 
the two-qubit state to encode one effective qubit: parallel 
or anti-parallel spins. The weak measurement technique 
described above will then partially collapse the state to- 
ward the parallel state under a null-measurement (no sin- 
gle electron tunneling). If now a 7r-pulse is applied to one 
of the spins with electron spin resonance, followed by a 
second null-measurement, this was shown to uncollapse 
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the state of the effective qubiti^ 

B. Optical polarization qubit 

Another interesting development is the experimental 
implementation of wavefunction uncoUapse for optical 
qubits using the polarization degree of freedom of single 
photons by Kim et The weak measurement was im- 
plemented by passing the photon through a glass plate 
oriented at the Brewster angle. Only the vertical po- 
larization is reflected off of the glass plate (with some 
probability). By placing a single-photon detector where 
the photon would have gone had it reflected, a null-click 
measurement partially collapses the polarization state to 
the horizontal polarization. The strength of the mea- 
surement can be increased by placing a series of plates 
in a row, effectively increasing the net probability of a 
vertically-polarized photon reflecting at some point. 

The wavefunction uncoUapse is done by inserting a 
half-wave plate (exchanging the amplitude of horizontal 
with vertical polarization) , and having the same number 
of plates traversed by the photon again. If none of the 
single-photon detectors click, the polarization state is un- 
coUapsed. This has been verified^ with quantum state 
tomography (with polarizer and single-photon counter 
placed after all of the reflecting plates) on the photon, 
conditioned on none of the other photon detectors fir- 
ing. The experiment showed an uncoUapsing fidelity of 
above 94% for measurement strengths up to 0.9. It was 
also pointed out that the information from the first weak 
measurement can be used for developing guessing strate- 
gies about the unknown initial state. Two such strate- 
gies were presented, and one was shown to be optimal. 
Of course, in the case where the measurement was subse- 
quently undone, these strategies did no better than ran- 
dom guessing. 

VII. CONCLUSION 

We have reviewed and extended recent developments in 
the theory (and experiment) of wavefunction uncoUapse 
by undoing quantum measurements. We have formulated 
the problem of wavefunction uncoUapse in terms of a con- 
test between the uncoUapse proponent, Plato, and the 
uncoUapse skeptic, Socrates, monitored by the arbiter, 
Aristotle. Plato claims to have the ability to uncoUapse 
wavef unctions, and this ability can be tested under the 
rules of the contest set forth. 



We have discussed several general features of the un- 
coUapse process in the abstract case, such as the upper 
bound on the success probability and quantum informa- 
tion aspects of the problem. In order to probabilistically 
undo the measurement, it is necessary to erase the in- 
formation extracted about the state in the first measure- 
ment. This is a necessary condition to uncoUapse the 
wavefunction, because otherwise various paradoxes arise. 
However, the information erasure is surely not a sufficient 
condition: the unitary evolution should also be properly 
reversed and, as the most experimentally challenging con- 
dition, the process should not bring decoherence, which 
requires a very good (ideal) detector. 

In addition to discussing the theory of wavefunction 
uncoUapse in the abstract case, we have also considered 
a variety of solid-state implementations and specific prac- 
tical strategies for wavefunction uncoUapse. Not only the 
success probability, these specific systems also allow the 
calculation of other characteristics of the process, such as 
the waiting time distribution for wavefunction uncoUapse 
in the charge qubit case. The cases of the charge qubit 
(with and without Hamiltonian dynamics), the phase 
qubit, and many entangled charge qubits have been ex- 
amined in detail. Additionally, we have also discussed 
two experimental realizations of this physics, based on 
the phase qubit and the polarization qubit, both of which 
have clearly demonstrated wavefunction uncoUapse with 
high fidelity. 

The ideality of the detector is necessary for perfect 
uncoUapsing, and we have only dealt with these kinds 
of detectors in the theory section of this paper (by this 
we mean the detector adds no extra decoherence to the 
system). If a detector is slightly non- ideal, then even 
a perfectly executed uncoUapse strategy will result in a 
slight infidelity in the final state. This is indeed the case 
in the experiments mentioned above although the fidelity 
was quite high. In such a situation there are two char- 
acteristics to contend with: the fidelity of uncoUapsing 
as well as the probability of claimed success. It is an 
open topic for future research how these characteristics 
are related. 
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